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Abstract 

Different formulas relying measurable fragment isotopic observables to the symmetry energy of 
excited nuclei have been proposed and applied to the analysis of heavy ion collision data in the 
recent literature. In this paper we examine the quality of the different expressions in the framework 
of the McGill Canonical Thermodynamic Model. We show that even in the idealized situation of 
canonical equilibrium and in the absence of secondary decay, these formulas do not give a precise 
reconstruction of the symmetry energy of the fragmenting source. However, both isotopic widths 
and isoscaling appear very well correlated to the physical symmetry energy. 

PACS numbers: 25.70Mn, 25.70Pq 
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I. INTRODUCTION 



One of the greatest challenge of intermediate energy heavy ion physics consists in the 
expeHmenta. measurement of the modifications to the sym„,etry energy c.,„ induced by 
density and temperature conditions different from the ones of ground state nuclei [H]. For the 
different astrophysical applications linked to the evolution and structure of compact stars 2], 
Q], the symmetry energy behavior at density far from saturation is of outmost importance. 
The high density behavior of the isovector equation of state is almost not constrained by 
experimental observations, and considerable uncertainties exist also in the behavior of the 
symmetry energy at subcritical densityj4]. In this regime, nuclear matter is unstable with 
respect to phase separation, mean-fie^ 



clusterization has to be considered 



d estimation can become severely incorrect [5| , and 
,|7|. Another important point concerns the temperature 
dependence of the symmetry energy, which is schematically treated or even neglected in the 
supernova explosion and proto-neutron star cooling modelizations 5|, ISj. 

A possible approach to this problem consists in comparing selected isospin observables 
to the output of a transport model where the isovector part of the equation of state can 
be varied 9| . This strategy has recently lead to very stringent constraints on the symme- 



try energy behavior 



101 ]. that a ppe ar reasonably consistent with the experimental results 



extracted from collective modes 



111 ] , nuclear massses 121] and neutron skins 131] . The draw- 



back of these analyses is that the results are model dependent, and different models do not 



produce fully consistent results 
finite temperature behavior. 

Alternatively, a simple formula has been proposed 



moreover no information can be extracted about the 



151 ] to extract directly the symmetry 



energy from experimental cluster properties obtained in the fragmentation of two systems 
of charge ^1,^2, mass Ai, at the same temperature T: 

.0 



T 



a 



{Zl/AD - [Zl/AD 



(1) 



where a is the so-called isoscaling parameter that can be measured from isotopic yields 16| . 
Applications of this formula to different fragmentation datajl^, 13] show that the measured 
symmetry energy tends to decrease with increasing collision violence. 

However eq. ([H) is not an exact expression. It was derived in the framework of macroscopic 
statistical models [l^, where many-body correlations are supposed to be entirely exhausted 



2 



by clusterisation, and it appears to be strongly affected by conservation laws and combina- 



torial effects 



19 



20| : secondary decays may strongly affect the value of a[2l|; finally, the c 



22 1, which is 



coefficient appearing in eq.([T]) should correspond to the symmetry free-energy 
equivalent to the symmetry energy only in the T ^ limit. 

In particular, S. Das Gupta et al. have shown that the difference in the neutron chemical 
potential A/i„ increases with the temperature if Csym is taken as a constant 19|]. In the grand- 
canonical ensemble, the difference of chemical potential between two sources at a common 
temperature T is linked to the isoscaling parameter by 



A/i„ = aT 



(2) 



This means that eq.([T]) cannot be exact at finite temperature. Experimentally the mea- 
sured value of aT decreases with increasing incident energy and/or collision violence. Let 
us suppose that the grand-canonical equality ([2]) is true also in the data, and that ncreasing 
collision violence does indeed correspond to increasing temperature. Then this implies that 
the physical symmetry energy coefficient explored in fragmentation data is not constant. 



but it is decreasing more steeply t 



with the interpretation o 
model calculations of ref. 



ref. 



17 



lan obtained in the rough data 17|, ll8| . This is consistent 
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23 



24 



251]. But it is inconsistent with the statistical 



26|, where a constant input symmetry energy coefficient produces 



an apparent Csym from eg. (HI) qualitatively coherent with the experiment. The statistical 



model MMM used in ref. 26| is similar, but not identical, to the statistical model CTM 
used in ref. jiol, w hich raises once again the question of the model dependence of the results. 
However, refs. [l9| , 26| do not compute the same observables either, and use different pre- 
scriptions for the symmetry energy. In these conditions, it is difficult to understand the real 
origin of the observed discrepancy. 

To progress on this issue, we present in this paper calculations with CTM made under 
similar conditions with respect to MMM. We will show that the two models produce very 
similar results. Moreover, we address the question of an "improved" formula which would 
be valid out of the T — > limit. We will show that, in the framework of this model, none of 
the different formulas proposed in the literature allows a reliable direct measurement of the 
symmetry energy. However, both the isoscaling observable and isotopic widths appear very 
well correlated with the physical symmetry energy, implying that ratios of these isotopic 
observables measured in different systems should allow to extract the physical trend. 
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II. THE MODEL 



The canonical thermodynamic McGill model is based on the analytic evaluation of the 
canonical partition function for the fragmenting source with A nucleons and Z protons 
(neutron number N = A — Z) aX a. given temperature T written as 



Q 



A,Z 



En 



n 



(3) 



where the sum is over all possible channels of break-up which satisfy the conservation laws; 



n. 



is the number of this composite in the given channel, and uja,j is the partition function 



of one composite with nucleon number a and proton number j: 



/i3 



(4) 



Here ma is the mass of the composite and Vf = V~Vq is the volume available for translational 
motion, where V is the volume to which the system has expanded at break up and Vq is 
the normal volume of A nucleons and Z protons. Concerning the choice of Zaj{int) used 
in this work, the proton and the neutron are fundamental building blocks thus Zi^Q{int) = 
zi^i{int) = 2 where 2 takes care of the spin degeneracy. For deuteron, triton, ^He and ^He 
we use Za,j{int) = {2saj + 1) exp{—i3eaj{gr)) where (3 = 1/T,ea,j{gr) is the ground state 
energy of the composite and {2saj + 1) is the experimental spin degeneracy of the ground 
state. Excited states for these very low mass nuclei are not included. For mass number 
a = 5 and greater we use the liquid-drop formula: 



^ajiint) = exp—[Woa - a(T)a 



2/3 



r 



T'' " ^ ' a^/^ "a eo 

The expression includes the volume energy, the temperature dependent surface energy, the 
Coulomb energy, the symmetry energy and contribution from excited states in the continuum 
since the composites are at a non-zero temperature. In this paper we will try different pre- 
scription for the symmetry energy coefficient, namely the same mass dependent prescription 
employed in the MMM model 26 1 



[a - 2jf T-a 

^sym ' I 



(5) 



-'sym 



CiC'}; 



CiCgd ^ , 



(6) 



with Ci = 1.7826, = 15.4941 MeV, Cg = 17.9439 MeV, or alternatively a more sophisticated 

271 ] , accounting for the vanishing of all surface 



surface and temperature dependent expression 



contributions at the critical point: 



■ rp2 _ rp2\ 5/4 



c 



sym 



a,T) = co-QT,a-i/M-^--- . (7) 



2^2 _|_ J^2 



where Cq = 28.165 MeV, and = 18 MeV is the critical temperature. To test the sensitivity 
of the different observables to the symmetry energy, a schematic constant coefficient will also 
be used. 

In using the thermodynamic model one needs to specify which composites are allowed in 
the channels. For mass numbers a=5 and 6, we include proton numbers 2 and 3 and for 
mass number a=7, we include proton numbers 2,3 and 4. For a > 8, we include all nuclei 
within drip-lines defined by the liquid-drop formula. 

The Coulomb interaction between different composites is included in the Wigner-Seitz 



approximation 27|, l28|]. For further details, see ref. 27 1. 



III. SYMMETRY ENERGY EVALUATIONS 

In this section we present the results of calculations made with the CTM model with 
the aim of reconstructing the input symmetry energy of the model from measurable cluster 
observables. When not explicitely stated, we will consider an excited fragmented source 
composed of Z = 75 protons with two different mass numbers A = 168, A = 186. This 



specific choice of source size was already employed in previous works 19|,|29|. Tabledgives the 



value of the isoscaling parameter obtained in the model and the resulting apparent symmetry 



energy from eq.([T]) for different values of the temperature and the break-up volume 
isoscaling parameter a is the value extracted from the slopes of differential cluster yields 



'he 



averaged over Z = 1,2, 3, 4, 5, similar to the procedure employed in the analysis of heavy ion 
data 9]. For these light isotopes, an excellent isoscaling is observed in the model 



3 

ion 

19(. The 



input symmetry energy in this exploratory calculation is fixed to Cgym = 23.5 MeV. 

The isoscaling parameter a decreases with increasing temperature independent of the 



break-up volume. This is in agreement with the results of previous works [19|, l26|| , as well 
as with the experimental observation of a decreasing a with increasing collision violence. 
Concerning the extracted values of Csym, an important dependence on the break-up volume 
is observed. Fot small break-up volumes, the apparent c^^m monotonically increases with 
temperature as already observed in a preceeding study with the CTM model where a was 
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Temperature 




a 


syTfi 


7.5 MeV 


4 


0.511 


26.10 MeV 


6.5 MeV 


4 


0.557 


24.65 MeF 


5.5 MeV 


4 


0.606 


22.72 MeV 


4.5 MeV 


4 


0.703 


21.56 MeF 


3.5 MeV 


4 


0.870 


20.75 MeF 


7.5 MeV 


6 


0.462 


23.6 MeV 


6.5 MeF 


6 


0.514 


22.75 MeF 


5.5 MeV 


6 


0.578 


21.65 MeV 


4.5 Mey 


6 


0.673 


20.63 MeF 


3.5 MeV 


6 


0.923 


22.02 MeF 



TABLE I: Isoscaling parameter averaged over Z = 1 — 5 and apparent symmetry energy from 
sq-dH) for a fragmenting source with Z = 75 at different temperatures and break-up volumes. The 
input symmetry energy in the model is Cgym = 23.5 MeV independent of the temperature. 

deduced from the chemical potential via the grancanonical expression eq.([2]) and not de- 
duced from the slope of fragment yields[19|. For higher volumes, the apparent Csym initially 



decreases as in ref. 
regime may not be explorec 
increasing collision violence 



261 ] and in the data, and it increases again when a saturates. This second 
in the data because the break-up temperatures saturate with 



171 ]. The same may be true for the analysis of ref. 26| which 



was made in the microcanonical ensemble; indeed in this latter the temperature does not 
increase linearly with excitation energy, producing a saturation in the apparent Cgym- The 
existence of these different behaviors shows that grancanonical formulas have to be handled 
with care: in the canonical or microcanonical models the slope of isotopic yields may not be 
directly related to the chemical potential. 

For the higher break-up volume, the input symmetry energy coefficient is well recovered 
at the lowest temperature. This was expected since eq.([T]) has been derived in the limit of 
vanishing temperature. Surprisingly, this does not seem the case if the break-up volume is 
small. In this case, the limit may be attained at lower temperatures where our fragmentation 
model cannot be safely applied any more. However, for all the situations considered, which 
cover a large range of thermodynamic conditions typically accessed in fragmentation exper- 
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iments, the deviation between the input Cgym and the approximation extracted from eq.([T]) 
never exceeds 12%, which is a reasonable precision considering the inevitable error bars 
induced by efficiency, event selection and thermometry in heavy ion collision experiments. 

To further progress on this analysis, we plot on the left part of FigH] the isoscaling a 
parameter as a function of the temperature with different choices of the symmetry energy 
parametrization. In all decreasing isoscaling parameter is found. The middle part 

of the same figure shows the resulting symmetry energy coefficient obtained by applying 
eq.([T]). We can see that the functional form of Cgym does not affect strongly the trend of the 
results. In particular, a decreasing a does not necessarily imply a decreasing physical sym- 
metry energy. Moreover, the temperature and surface dependence of the physical symmetry 
energy affects the predictive power of eq.([T]). In no case the extracted coefficient approaches 
the symmetry energy of the fragments used for the isoscaling analysis, however at a given 

the input symmetry energy of the 



SOj. This means that the isoscaling 



value of temperature, it qualitatively follows the trend o 
fragmenting source, as expected in the Weisskopf regimejc 
properties of the lightest fragments appear well correlated to the symmetry energy of their 
emitting source, even out of the evaporation regime. 

From the observations of Fig{T]and tabled we can already draw some partial conclusions. 
An important point concerns the fact that observing a or c^y^ decreasing with the collision 
violence cannot be taken as an evidence that the physical symmetry energy does so. Only 
a detailed comparison with a model may allow to extract the physical symmetry energy. 
Different models have to be very carefully compared to a large set of independent observables 
before one can extract any conclusion. As a second remark, both MMM and CTM models 
tend to agree that at low temperature eq.([T]) gives a good reproduction of the physical 
symmetry energy. 

This means that results obtained from intermediate impact parameter collisions in the 
neck region 3ll,l32] (where in principle the matter is at low density but also relatively cold), or 
analyses of quasi-projectiles produced in peripheral collisions 17|] (where the system is close 
to normal density and the temperature behavior could be disentangled from the density 
behavior) are better suited to this study than central collisions in the multifragmentation 
regime. 

To progress on the issue of the determination of the in-medium modifications to the 
symmetry energy, it would be extremely useful if we could have a formula more adapted 
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to the finite temperature case. To this purpose, we turn to check two other expressions 
proposed in the literature to access the symmetry energy from fragment observables. 



IV. THE INFLUENCE OF FRACTIONATION 



In nuclear multifragmentation reactions, the asymmetry term influences the neutron- 
proton composition of the break-up fragments. 

Interpreting multifragmentation in the light of first-order phase transitions in multicom- 
ponent systems, the neutron enrichment of the gas phase with respect to the liquid phase 
comes out as a natural consequence of Gibbs equilibrium criteria and a connection between 
phases chemical composition and the symmetry term can be established 33|, [3^. Interest- 



ing enough, the phenomenon of isospin fractionation which is systematically observed in 



analyses of multifragmentation data [15 



35 



36, 



37|, seems to be a generic feature of 



phase separation independent of the equilibrium Gibbs construction [38|. Indeed, dynamical 



models of heavy ion collisions 



, y] where fragment formation is essentially ruled by the out 
of equilibrium process of spinodal decomposition also exhibit fractionation. Adopting an 
equilibrium scenario for the break-up stage of a multifragmenting system, Ono et al. 39 1 
derive an approximate grandcanonical expression which connects the symmetry term with 
the isotopic composition of fragments obtained in the break-up stage of two sources with 
similar sizes in identical thermodynamical states and differing in their isospin content. 



(8) 



under the hypothesis that the isotopic distributions are essentially Gaussian and that the 
free energies contain only bulk terms. Here, a{j) is the isoscaling slope parameter of a 
fragment of charge j, and < a >j stands for the average mass number of a fragment of 
charge j produced by the source i{= 1, 2) at the temperature T. 

In the limit of vanishing temperature, fractionation can be neglected and j/ < a can 
be replaced by the corresponding quantity of the sources Zi/Ai ISj giving back eq.([I]). In 
the opposite case, fragment yields are predicted to be sensitive to their proper symmetry 
energy and not to the symmetry energy of the emitting source. 

Figure [2] gives, as a function of the cluster atomic number j, the apparent symmetry 
energy coefficient extracted from eq.([8]) for different conditions of temperature, free volume. 
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source isotopic content, and source size. In all cases the input symmetry energy has been 
fixed to the constant value Csym = 23.5 MeV. We can see that eq.® leads to a global 
systematic overestimation of the input symmetry energy. The response of the different 
fragments depends on their size. Non realistic values are obtained for the lightest fragments. 
The lightest fragments isotopic distribution j/ < a > is very sensitive to the number of 
isotopes considered in the calculation, which show strong binding energy fluctuations. It is 
not surprising that they cannot be treated by eq.®, which implies a behavior dominated 
by the bulk for all fragments. Turning to the increase in the apparent symmetry energy 
for the heaviest fragments, this is most probably due to the failing of the grancanonical 
approximation in eq.([H]) when the fragment size becomes comparable to the source size, as 
previously discussed in ref.j2^. One should also note that for heavy fragments isoscaling 
tends to be violated in the model[l^, and the determination of an isoscaling slope becomes 
largely arbitrary. Clusters of charge j > 5 and smaller than approximately one tenth of the 
source size are best suited to this analysis. If we limit ourselves to such intermediate mass 
fragments, we can see that the apparent Cgym coefficient is reasonably independent of the 
available volume, source isospin and mass. A temperature dependence is still apparent and, 
as in the case of FiglH does not exceed 10%. These results are in good agreement with the 
findings of ref. 26| in the framework of the MMM model. 

Figure [3] shows the response of eq.® to a symmetry energy depending on the temperature 
and on the fragment size through eq.(l7l). The behavior is very similar to the one displayed 
in Figure [2] above. This means that it is not possible to extract the surface dependence 
by looking at the behavior as a function of the charge. The temperature dependence for 
a fixed charge (right part of Figure [3]) conversely shows a good correlation, meaning that 
the temperature dependence could be extracted studying the isoscaling for a given charge 
at different excitation energies. 



An alternative expression has been derived in ref. [40|| connecting the symmetry energy of 
a cluster of size a to the width of its isotopic distribution. Indeed, a Gaussian approximation 
on the grandcanonical expression for cluster yields gives 

where cr|(a) indicates the width of the isotopic distribution of a cluster of size a, and / = 
a — 2j. In principle the Csym coeffcients appearing in eqs.® and ([9]) correspond to symmetry 



free-energy coefficients, that is they include an entropic contribution. However if we neglect 
the / dependence of the excitation energy and entropy associated to a given mass a, they 
coincide with Cgym. defined by eq.dS]). 

Figini shows the apparent Cgym extracted from the fluctuation formula eq.([n]) as a function 
of the cluster charge j at different temperatures and isospin values. This observable shows 
a linearly increasing behavior similar to the one displayed by eq.®, but the overall quality 
of reproduction is improved. 

In this picture the input symmetry energy was taken as a constant. To see if eqs.QjQ 
can be used to extract from observable cluster data the possible surface/temperature/density 
dependence of the symmetry energy, we additionally show in Figure \5[ for a cluster charge 
j = 10, the apparent Cgym extracted from eqs.dH]),® as a function of the input Csym- We can 
see that both formulas produce a bias. This implies that a quantitative estimation of the 
symmetry energy coefficient cannot be obtained from these expressions. We recall that the 
presented calculation completely neglects secondary decay, which is expected to considerably 
increase this bias. 

If the quantitative values do not match, we observe however a good linear correlation in 
both cases. This means that the analyzed observables show an excellent sensitivity to the 
isovector equation of state. In particular, the relative variation of the extracted symmetry 
coefficient should be reliable, and it would be very important to check whether the linear 
dependence survives to secondary decay. We recall that in the framework of the MMM 
model which (at variance with CTM) contains an afterburner, the proportionality is kept 
also after secondary decay [4^. 

V. CONCLUSIONS 

In this paper we have studied the sensitivity to the symmetry energy of different isotopic 
observables measurable in heavy ion collisions, in the framework of the McGill Canonical 
Thermodynamic Model. 

We conclude that, even in the idealized limit of thermal equilibrium, no isotopic ob- 
servable can allow to reconstruct the physical symmetry energy of the excited fragmenting 
source in a model indpendent way. The different models have to be very carefully compared 
to a large set of independent observables before one can extract any conclusion. In the low 
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temperature limit, the well spread expression eq.([T|) gives a good measurement of the sym- 
metry energy of the source, and is never sensitive to the symmetry energy of the fragments. 
At high temperature in the multifragmentation regime, no formula gives a satisfactorily 
reproduction of the input Csy^- 

However we confirm, in agreement with the results of previous studies with different 



models 



17 



19 



26 



40| that both the isoscaling variable and isotopic widths show a 



very strong sensitivity to the strength of the symmetry energy. 

This means that the use of the different formulas proposed in the literature should allow to 
better constrain Csym- In particular, if a drastic reduction of the effective Cgym in the nuclear 
medium is observed, we should be able to see it by calculating the differential response of 
eqs. ([H]) , © for different excitation energies. 

It is important to stress that this study neglects secondary decay which may have dramatic 
consequences on isotopic observable [2 ll|. It is clear that the issue should be investigated. 



Moreover differential observables 
should be analysed. 



iblej2]|. 



where the effect of secondary decay may cancel out 
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FIG. 1: As a function of the temperature, the left side shows the isoscaUng a parameter, the 
middle part shows the resulting apparent Csym from eq.([T]), and the right part the ratio between 
the apparent Csym and the input Csym of the source. Upper lines (black online): constant Cgym = 23.5 
MeV. Middle lines (red online): mass and temperature dependent Cgym. from eq.([7|). Lower lines 
(green online): mass dependent Cgym from eq.®. The lines with symbols in the middle panel give 
the input symmetry energy for a mass equal to the source(s) mass. The freeze-out volume is fixed 
to V/Vo = 6 
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FIG. 2: Apparent symmetry energy calculated from eq.® as a function of the cluster charge. 
Upper left: calculations at different temperatures T = 4 MeV (dotted line), T = 5 MeV (full line), 
T = 6 MeV (dashed line) for a source charge Z = 75, source masses Ai = 168, A2 = 186, and a 
volume V/Vo = 6. Upper right: calculations at different volumes V/Vq = 3 (dashed line), V/Vq = 4 
(dotted line), V/Vq = 6 (full line) for a source charge Z = 75, source masses Ai = 168, A2 = 186, 
and a temperature T = 5 MeV. Lower right: calculations at different isospin values Ai = 177 
and A2 = 168 (dotted line), Ai = 186 and A2 = 168(full hue), Ai = 186 and A2 = 177 (dashed 
line) for a source charge Z = 75, a temperature T = 5 MeV, and a volume V/Vq = 6. Lower left: 
calculations for different source sizes Z = 75, Ai = 186 and A2 = 168 (full line), Z = 60, Ai = 149 
and A2 = 135(dashed line), Z = 90, Ai = 222 and A2 = 202 (dotted line) for a source temperature 
T = 5 MeV, and a volume V/Vq = 6. 
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FIG. 3: Left side: Apparent symmetry energy calculated from eq.([8]) as a function of the cluster 
charge for a volume V/Vq = 6. Left side: Temperature T = 4 MeV and the input Csym displayed 
by dotted line is taken from eq.®. Right side: same as left, but for a fragment charge Z = 10 as 
a function of the temperature and the input Csym displayed by dotted line is taken from eq.©. 
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FIG. 4: Left side: apparent Csym extracted from the fluctuation formula eq.(l9|) as a function of the 
cluster charge j for a source charge Z = 75, source mass Ai = 168, and a volume V/Vq = 6, and 
temperature T = 4 MeV (full line), T = 6 MeV (dotted hne), T = 8 MeV (dashed line). Right 
side: same as left, but for a temperature T = 6 MeV and source mass Ai = 168 (full line) and 
A = 186 (dotted hne). The x axis is j/A. 
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FIG. 5: apparent Csym extracted from eq.® (left side) and eq.([8]) (right side) as a function of the 
input Csym for a cluster charge j = 10. Conditions are as in Figs l2|4l above. 
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